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THERMTRAJ: A FORTRAN PROGRAM TO COMPUTE THE 
TRAJECTORY AND GAS AND FILM TEMPERATURES 
OF ZERO PRESSURE BALLOONS 

By 

Walter J. Horn 
and 

Lei and A. Carlson 
Texas A&M University 

SUMMARY 

A FORTRAN computer program called THERMTRAJ is presented which 
can be used to compute the trajectory of high altitude scientific zero 
pressure balloons from launch through all subsequent phases of the 
balloon flight. In addition, balloon gas and film temperatures can be 
computed at every point of the flight. The program has the ability to 
account for ballasting, changes in cloud cover, variable atmospheric 
temperature profiles, and both unconditional valving and scheduled, valv- 
ing of the balloon gas. The program has been verified for an extensive 
range of balloon sizes (from 0.5-41.47 million cubic feet). Instructions 
on program usage, listing of the program source deck, input data deck and 
printed and plotted output for a verification case are included. 

INTRODUCTION 

The program described in this report was developed under NASA 
Contract NAS6-3072 to accurately compute the trajectory of high altitude 
scientific zero pressure balloons, and to compute the balloon film and 
gas temperature. A set of five ordinary differential equations based 



upon the dynamics and thermodynamics of the balloon system have been 
formulated. These equations are solved using a fourth order Runge- 
Kutta solution scheme with a variable time step. The program allows 
the user to specify a ballast schedule, to select from five atmospheric 
models (or supply his own) to adjust for changes in cloud cover and 
to account for balloon gas valving. 

The background analysis for the computer program is contained in 
the reference, "A Unified Thermal and Trajectory Model for the Prediction of 
High Altitude Balloon Performance." The purpose of this report is to 
detail the use of the computer program. The following section is a 
presentation of the governing ordinary differential equations and the 
details of the options available to the user of the program. Subsequent 
sections will deal with the input data requirements and formats, the 
output format (both printed and plotted), a program listing, and the results 
of a verification study conducted during the development of the program. 

PROBLEM DESCRIPTION 

The trajectory of a balloon system is depends upon the dynamics of 
the system and the thermal contributions from its atmospheric environment. 
Figure 1 contains a schematic of a typical zero-pressure balloon system 
indicating the source of thermal loading on the system. The results of 
the analyses are reported in the Reference. They show that the vertical mo- 
tion of the balloon is governed by the following set of five ordinary 
differential equations: 


z , dz/dt 
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Symbols: T = Energy transmitted 

R = Energy reflected 
A = Energy absorbed 
m = Mass of payload plus ballast 

r 

nif - Mass of balloon film 
m g = Mass of balloon gas 


Figure 1. Schematic of typical balloon system 
with the thermal inputs to the balloon. 
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where; t = time 

v = vertical velocity of the balloon system 
g = acceleration due to gravity 
p n = density of air 

a 

V g = volume of the balloon 

iry = mass of the balloon film 

m ~ mass of gas 
y 

m = mass of the payload = mass of the balloon system - (m f + m ) 
h * y 

C Q = coefficient of drag 

A = effective balloon cross-sectional area = 1.2089 V 

y 

z = balloon altitude 
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Pg = balloon gas pressure 


OFpoq}J* A ®Z IS 
Mg = molecular weight of balloon gas QUALITY 

R = universal gas constant 

Tg - temperature of the balloon gas 

§ g = volume flow rate of gas exhausted at entry into float 
% = mass flow rate of balloon gas during controlled valving operation 


T f = temperature of balloon film 
Q.p = net heat flux to film 
Cf * specific heat of balloon film 
C = specific heat of balloon gas 

H g 

Qq = net heat flux to balloon gas 


M a = molecular weight of air 
T g = temperature of the air 


and the radiative heat transfer to the balloon film can be expressed as 

q F =0 /4 G V eff + G Vw.eff + -int CT C T g 4 ' V> + 


CHgf < T g - V + CH fa< T a ‘ V " Wf - W> 


( 6 ) 


while the radiative heat transfer to the balloon gas can be expressed as 
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Where the effective coefficient of absorptivity of the balloon film, the 
effective coefficient of infrared interchange between the balloon film 
and gas, the effective coefficient of emissivity of the balloon film, the 
effective coefficient of absorptivity of the balloon gas, and the effective 
coefficient of emissivity of the balloon gas are expressed in terms of 
the gas and film radiative properties as: 

T w,sol (l-a q ) 

1 - r.. _ fl-aj 


a,, ~ a,. 

w,eff W 


1 + 
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where G = solar constant 

r g = earth reflectivity (albedo) 

a = Stefan-Bol tzman constant 

CH f = convective heat transfer coefficient between the balloon film 
gT and gas 

CH-, 25 convective heat transfer coefficient between the balloon film 
Ta and air 



= black ball temperature 


S = balloon surface area = 4.835976 V ^ 

a = coefficient of absorptivity of the balloon film in the infrared 
w spectrum 
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t.. cat s coefficient of transmissivity of the balloon film in the solar 
w,SOi spectrum 

a = coefficient of absorptivity of the balloon gas 

r cm- coefficient of reflectivity of the balloon film in the solar 
w,so “ spectrum 

e = coefficient of emissiv^y of the balloon gas 
y 

e = coefficient of emissivity of the balloon film 
w 

r = coefficient of reflectivity of the balloon film in the infrared 
w spectrum 

t = coefficient of transmissivity of the balloon film in the infrared 
w spectrum 

Differential equations 1 through 5, along with the defining equations 
6 through 12, are solved using a Kutta-Simpson, one-third rule, fourth order 
Runge-Kutta solution technique. 

The atmospheric properties are computed using any of five "standard" 
atmosphere models or a user supplied model. The five standard models are 
the 1962 Standard Atmosphere model and four models based upon NSBF temper- 
ature-altitude profiles measured during the winter, spring, summer and fall 
of the year. 

Ballasting is accomplished by providing a ballast rate-time schedule 
corresponding to the actual ballast schedule of the balloon flight being 
simulated. This schedule represents the ballasting rate at all times 
during the flight; therefore, each ballasting operation is represented with 
an entry at both the beginning and the end of the ballast drop. 

Controlled balloon gas valving operations are modelled in a manner 
similar to that for ballast drops. That is, a valving mass flow rate-time 
schedule is provided for every point of the flight. Therefore, every val- 
ving operation is represented by two entries in the valving schedule; one 


7 



at the beginning of the valving operation and one at the end. Valving of 
the balloon gas at the entry into float is handled automatically by the 
program by monitoring the volume of the balloon. When the balloon volume 
exceeds the maximum volume of the balloon during a time step in the solution, 
a sufficient volume of gas is expelled to prevent the volume of the balloon 
from exceeding the maximum volume of the balloon and a smaller time step is 
used in the solution. 

Cloud cover conditions are specified by a time schedule of albedo values 
and black-ball temperature prefile adjustments at every point in the flight. 
Therefore, if accurate cloud cover information is available for the balloon 
flight being simulated it can be modelled within the program. 

A variable time step solution of the governing differential equations 
was developed to economize the solution scheme. Basically, the procedure fol- 
lowed was to double the time step if the computed values of the increment of 
the balloon altitude, velocity, and film temperature and gas temperature were 
less than some user provided tolerance for the previous three steps in the 
solution. In a similar manner, the time step was reduced by a factor of 0.5 
if the computed values of the same increments mentioned above exceed the user 
provided tolerance during any step of the solution and that solution step was 
repeated. 

PROGRAM USAGE 

The computer program is written in FORTRAN IV programming language and 
is operational on the Texas A&M Amdahl 470 V/6 and Amdahl 470 V/7B computers. 

The program uses plotting subroutines developed by California Computer Products, 
Inc. (Cal -Comp) which are available on the Texas A&M computer system. 
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This section gives an outline of the procedure for using the 
THERMTRAJ program *or solving the equations, developed in Reference 1 
and summarized in the previous section of this report, governing the 
vertical motion of the balloon system. Table 1 contains a summary of 
the input data cards necessary to run a flight simulation and Table 2 
contains a summary of the units of the input data. Table 3 contains a 
list of the default values of the input variables for the program. Only 
those variables that are to be modified need to be assigned values on 
the data card associated with the READ (5, FIPUT) statement of line 91 

of the source program listing of Appendix A. Table 4 contains a list of 

the FORTRAN variable names used in the program and their definitions. 

Ta M o K r nrtfa ■? a pjimmaww n ■f f Ua aX 1 4*1%^ MMAAutM 

iMviw wr WUiivw iiio u auiidiiut jr vi uic v/uupul Ul l/lic ptuyioili. 

Figure 2 contains a typical four segment atmospheric temperature- 
altitude profile used to define the atmospheric properties at each point 
in the balloon flight and Figure 3 contains a summary plot of the five 
optional profiles that are contained in the THERMTRAJ program. A typical 
three segment black-ball temperature-altitude profile is presented in Figure 
4 indicating the FORTRAN symbolic names of the coordinates of each of the 
break-points on the profile. 

The user may begin the balloon simulation at any point of the balloon 
flight. The following input data in the NAMELIST FIPUT must be modified to 
correspond to the conditions which exist at the point the simulation is 
initialized: 

TIMOL - Time of initialization of simulation (hrs) 

CUTDT - Length of flight from the point the simulation was initiated 
(hrs) 
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AS CRT and HFLIT - Must be adjusted so that 

HFLIT - BALT 

a§crt — 

yields either the time from initial start to the estimated 
time of float ojr the time when TBBC(I) is to be used instead 
of the three segment TBB profile 
MGAS - Mass of gas at the start of the simulation (g) 

MPAYI- Mass of payload, including remaining ballast at the start of 
the simulation (g) 

All other NAMELIST variables are as if the computations were started at the 
original launch time and altitude. 

The variables of the ballast drop schedule, the cloud cover schedule 
and the valving schedule must also be modified to correspond to the initial 
s*art time of the simulation. Specifically, the times in the arrays TBDRP(I), 
TCLDC(I), and TVALVE(I) must be times measured from the time of the initial- 
ization of the flight simulation, not from the launch time of the balloon. 

All input experimental data (EMB, EAIRT, EFILMT, EGAST) must be for 
times which occur after the time of the start of the simulation. 

Finally, the value of four variables within the body of the main pro- 
gram must be altered to correspond to the actual beginning of the simulation. 
The initial values of TR, TS, Zl, and Z2 (see lines 281-284 of the source 
deck of Appendix A) must be balloon gas temperature, film temperature, al- 
titude, and velocity at the time of the beginning of the simulation. 
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Figure 2 . Typical Atmospheric Temperature Model 




Figure 3. Five optional atmospheric temperature profiles available in the THERMTRAO program 
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Figure 4. Typical black-ball temperature profile 


Table 1. Summary of Input Cards 
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Read 

Order 

Variables 

Format 

1 

NTITLE 

20A4 

2 

NAMEL I ST /F IPUT/TIMOL , CUTDT ,BALT , ASCRT , HFL IT, 

Namel ist 


PASL,MASSF,MGAS,MPAYI,VOLMAX,DAYL,TWILT,ALBED, 



FLSOD , PT INC , DTMAX , SPHTW , SPHTG , MOL A , MOLG , ALFRES , 



EG, ALFSO.TWSOL , RWSOL , EW ,RW,TAUW,TASL , HTROP ,TATR0> 



HTR0P2 ,TATR02, HTROP3 ,TATR03,HFL0T,TAFL0,HBB2 ,HBB3, 



T2T,T3T,T4T,DPLT, IATMS 


3 

NBDRP 


4 

TBDRP(I) and BRATE(I) ; 1=1, NBDRP 
(omitted if NBDRP=0) 

2F10.2 

5 

NCLDC 

15 

6 

TCLDC(I) ,ALBC(I) and TBBC(I) 
(omitted if NCLDC=0) 

3F10.2 

7 

NVALVE 


8 

TVALVE(I) and VALVRT(I) 

2F10.2 


(omitted if NVALV'£=0) 


9 

NEX 

13 

10 

ELN1,ELN2,ELN3 

3F3.0 


(omitted if NEX=0) 


11 

TT(J),J=l,3j EMB(I) ,EAIRT(I) ,EFILMT(I) ,EGAST(I) 

3F3.0.1X, 


(omitted if NEX=0) 

4F10.4 
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Table 2 Summary of Units of Input Data 


TIMOL - 
CUTUT - 

hrs. 

hrs. 

ORIGINAL PAGE ISI 
OF POOR QUALITY 

HTR0P3 

TATR0P3 

- 

m 

°K 

BALT - 

m 


HFLOT 

- 

m 

ASCRT - 

m/ml n . 


TAFLO 


°K 

HFLIT - 

ni 


HBB2 

•u 

m 

PASL - 

millibars 


HBB3 

** 

m 

MASSF - 

g 


T2T 

- 

°K 

MGAS - 

g 


T3T 

m 

°K 

MPAYI - 



T4T 

- 

°K 

VOLMAX - 

m3 


DPLT 

m 

dimensionless 

DAYL - 

hrs. 


IATMS 

- 

dimensionless 

TWILT - 

hrs. 


NBDRP 

- 

dimensionless 

ALBED - 

dimensionless 


TBDRP 

- 

min. 

FLSOD - 

cal/mvmin. 


BRATE 

- 

g/min. 

PTINC - 

hrs. 


NCLDC 

- 

dimensionless 

DTMAX - 

min. 


TCLDC 


min. 

SPHTW - 

cal/g/°K 


ALBC 

«■« 

dimensionless 

SPHTG - 

cal/g/°K 


TBBC 



MOLA - 

dimensionless 


NVALVE 

- 

dimensionless 

MOLG - 

dimensionless 


TVALVE 

- 

min. 

ALFRES - 

dimensionless 


VALVki 

9m 

g/sec. 

EG 

dimensionless 


NEX 

- 

dimensionless 

ALFSO - 

dimensionless 


ELN1 


hrs. 

TWSOL - 

dimensionless 


ELN2 

- 

min. 

RWSOL - 

dimensionless 


ELN3 


sec. 

EW 

dimensionless 


TT(1) 

- 

hrs. 

RW 

dimensionless 


TT 2) 

to* 

min. 

TAUW - 

dimensionless 


TT(3) 

m 

sec. 

TASL - 

°K 


EMB 

- 

millibars 

HTROP - 

m 


EAIRT 

- 

°C 

TATRO - 

°K 


EFILMT 

- 

°C 

HTROP 2 - 

m 


EGAST 

*• 

°C 

TATROP2 - 

°K 
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Table 3. Summary of Input Default Values 




:: t, 

O' 


1=4 


pOU'^ 




TIMOL = 18.0 
CUTDT =24.0 
BALT = 120.0 
A SORT = 304.8 
HFLIT = 36700.0 
PASL = 1013.3 
MASSF = 1.842E5 
MGAS = 6.9221E4 
MPAYI = 1.9682E5 
V0LMAX = 6.6375E4 
DAYL = 12.0 
TWILT = 0.75 
ALBED = 0.18 
FLS0D = 1.9892E4 
PTINC = 0.01666 
DTMAX = 5.0 
SPHTW =0.55 
SPHTG = 1.24119 
MOLA = 28.9644 
MOLG = 4.0026 
ALFRES = 0.0026 
EG = 0.000312 
ALFSO = 0.001 


TWSOL = 0.885 
RWSOL =0.114 
EW = 0.031 
RW = 0.127 
TAUW =0.842 
TASL = 288.15 
HTROP = 11000.0 
TATROP = 215.65 
HTR0P2 =20000.0 
TATR02 = 216.65 
HTR0P3 = 32000.0 
TATR03 = 228.65 
HFLOT = 47000.0 
TAFLO = 270.65 
HBB2 = 11000.0 
HBB3 = 20000.0 
T2T = 214.4 
T3T = 214.4 
T4T = 214.4 
DPLT = -1.0 
IATMS = 5 
RGAS = 8.31432E7 
EDOTG =0.0 
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Table 4. List of FORTRAN variable names 
used within the program 
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FORTRAN 

Variable 

ABAL 

AG 

AGEF 

ALBC(I) - 
ALBED 
ALFRES - 
ALFSO 
A SORT 
ATM0S2 - 

AWEF 

AXIS1 

BALT 

BRATE(I) - 
C 

CD 

CFH 

CG 

CH 


Description 

2 

Surface area of balloon (m ) 

Acceleration due to gravity (m/mln ) 

Effective infrared absorptance of balloon gas 

Array of Albedo factors corresponding to times TCLDC(I) 

Current Albedo factor (0.18 for clear skies; 0.57 maximum) 
Absorptance of balloon gas 
Absorptance of balloon film 
Estimated balloon ascent rate (m/min) 

Subroutine to compute balloon altitude in meters when 
the altitude is given in millibars of pressure 
Effective infrared absorptance of balloon film 
Plotting subroutine (plots axes)(Local) 

Launch altitude (m) 

Array of ballast rates corresponding to times TBDRP(I) (g/min. ) 
Thermal conductivity of air (calorie/min/m/°K ) 

Coefficient of drag 

Function subprogram to compute the heat transfer coefficient 
between the balloon film and air 
Thermal conductivity of helium 

Heat transfer coefficient between balloon film and air 
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FORTRAN 

Variable 

CHF 

CHN 

CHR 

CHRES 

CIT 

CIITDT 

DAYL 

DINCH 

DPLT 

DTIME 

DTIMEN 

DTMAX 

DTMIN 

DX1 

EAIRT(I) 

EALT(I) 

EDOTG 

EDOTV(I) 


n . . ORIGINAL PAGE (S 

Description q F P0 OR QUALITY 

- Forced convection coefficient between balloon film and air 

- Natural convection coefficient between the balloon film and air 

- Heat transfer coefficient between the balloon film and gas 

- Function subprogram to compute the heat transfer coefficient 

between the balloon film and gas 

- Parameter used to determine relationship of the sun to the 

balloon (day/night) 

- Length of flight to be simulated (hrs.) 

- Length of daylight (hrs. input; converted to min. within program) 

- Time of day (hrs.) 

- Positive value of DPLT generates plotted output 

- Time increment (min.) 

- Minimum time increment (min.) 

- Maximum time increment (min.) 

- Plotting time scale increment (min. /inch of plot) 

- Plotting time scale increment (min. /inch of plot) 

- Array of measured ambient air temperature (°C) corresponding 

to the array of times ETIME(I) 

- Array of measured altitudes (m) corresponding to the array 

of times ETIME(I) 

- Gas volume flow rate for "burping" of gas when float altitude 

3 

is exceeded (m /sec) 

- Array of gas mass flow rate during valving operation (g/sec.) 

corresponding to the array of time TVALVE(I) 
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FORTRAN 

Variable 


Description 


« 
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EFilMT(I) - Array of measured balloon film temperature (°C) corresponding 
to the array of times ETIME(I) 

EG - Emissivity of the balloon gas 

EGAST(I) - Array of measured balloon gas temperature (°C) corresponding 

to the array of times ETIME(I) 

EGEF - Effective emissivity of the balloon gas 
EINT - Effective interchange emissivity between the balloon gas 
and film 

ELANCH - Launch time (sec.) = 3600.*ELNI+6Q.*ELN2+ELN3 

ELN1 - See ELANCH (hrs.) 

ELN2 - See ELANCH (min.) 

ELN3 - See ELANCH (sec.) 

EMB(I) - Array of measured altitude (millibars) corresponding to the 
array of times ETIME(I) 

EMB 1(1) - Array of In (pressure at launch si te/EMB ( I ) ) j plotting parameter 

ETIME(I) - Array of times for which measured temperature and altitude 

data have been input (sec) s 36Q0.*TT(1)+60.*TT(2)+T(3)- 
ELAUNCH 

EW - Emissivity of the balloon film 

EWEF - Effective emissivity of the balloon film 

EXPO - Exponent used in computing the pressure, density and temper- 

ature as a function of balloon altitude 
EXPT - Exponent used in computing the pressure, density and temper- 
ature as a function of balloon altitude 
EXP2 - Exponent used in computing the pressure, density and temper- 

ature as a function of balloon altitude 
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Variable 
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EXP3 - Exponent used in computing the pressure, density and temper- 
ature as a function of balloon altitude 
FK1 ,FK2,FK3,FK4,FL1 ,FL2,FL3,FL4,FM1 ,FM2,FM3,FM4,FN1 ,FN2,FN3,FN4,F01,F02, 
F03.F04 - Runge Kutta parameters of the general form: 


y i (t+dt)=y i (t)+l/6 (Fil+2Fi2+2Fi3+Fi4) 
for i = K,L,M,N,0 

FLIFT - Free lift on the balloon (kg) 

4 o 

FLSOD - Solar radiation constant (1.9892 10 cal/nr/min) 

FLXSD - Solar radiation flux adjusted for the time of day 

FMTOT - Total mass of balloon system minus the mass of gas (kg) 



F2 


F3 

F4 


> 


- Right hand sides of the five governing first order ordinary 
differential equations. 



GR - Grashof Number 

GY - Acceleration due to gravity (m/sec2) in Subroutine ATM0S2 

H - Parameter in subroutine ATM0S2 to designate altitude (m) 
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FORTRAN 

Variable 

HB 

HBAL - 
HBB2 - 

HBB3 - 

HFLIT - 

HFLOT - 
HMGAS - 

HTR 

HTROP - 
HTR0P2 - 
HTROP3 - 
HTS 

HVOLG - 


ORIGINAL page 

Description OF POOR QUALITY 

Parameter in subroutine ATM0S2 to designate the base altitude 
of the various layers of the standard atmosphere (m) 

Altitude of balloon (m) 

Black ball altitude corresponding to black ball temperature T2T 
in a black ball temperature-altitude three segment profile (m) 
Black ball altitude corresponding to black ball temperature T3T 
in a black ball temperature-altitude three segment profile (m) 
Estimated float altitude (m) used to determine TIMFL. When time 
exceeds TIMFL, TBBC(I) used instead of three segment TBB profile. 
Nominal float altitude (m) 

Duplicate storage location for the mass of gas during each time 
increment (g) 

Duplicate storage location for the balloon gas temperature during 
each time increment (°K) 

Altitude corresponding to the atmosphere temperature TATRO in a 
four segment temperature-altitude profile (m) 

Altitude corresponding to the atmosphere temperature TATR02 in a 
four segment temperature-altitude profile (m) 

Altitude corresponding to the atmosphere temperature TATR03 in a 
four segment temperature-altitude profile (m) 

Duplicate storage location for the balloon film temperature 
during each time increment (°K) 

Duplicate storage location for the balloon volume during each 
time increment (m 3 ) 
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FORTRAN 

Variable 

HZ1 

HZ2 

IATMS 


IPC 

ITC 

I VOL 

IVT 

KUTTA 

LINE 

LP 

MA 

MASSF 


Description 

Balloon altitude during each time increment (m) 

Balloon velocity during each time increment (m/sec) 

Parameter to choose atmospheric model, i.e., 

IATMS = 1 NSBF Winter Atmosphere 

2 NSBF Spring Atmosphere 

3 NSBF Summer Atmosphere 

4 NSBF Fall Atmosphere 

5 1962 Standard Atmosphere 

6 User supplied atmosphere 
Printer control character 

Parameter used in varying the ting ’ncrement of the 
solution of the governing differential equations 
Parameter used to identify regions of the flight that 
"burping" of balloon gas is required 
Parameter used to adjust the volume flow rate of the 
gas during the "burping" process 
Parameter used to direct the Runge-Kutta solution of 
the governing ordinary differential equations 
Local system subroutine used in the plot portion of the program 
Lapse rate of the layers of the 1962 standard atmosphere 
as contained in subroutine ATM0S2 
Molecular weight of air (Subroutine ATM0S2 parameter) 

Mass of balloon film (g) 
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FORTRAN 

Variable 


Description 


ORIGINAL PAG2 IS 
OF POOR QUALITY 


MGAS 

MGASD 

MOLA 
MOLG 
MPAY 
MPAY I 

MPD(I) - 

NBDRP 

NCLDC 


NEX 

NEX1 

NEX2 

NPC 

NPSUN - 
NSUN 

NTITLE(I) - 
NVALVE - 

NYC 

NYCP1 - 
NYCP2 - 
PA 


Mass of gas(g) 

Increment of the mass of gas for each time step of the 
solution 

Molecular weight of air (main program) 

Molecular weight of balloon gas 
Mass of the balloon payload (g) 

Mass of the balloon payload at launch including ballast (g) 
Array of mass of balloon payload associated with time of 
ballast drops TBDRP(I) 

(Number of ballast drops + 1) * 2; input number of ballast 
drops required 

(Number of cloud cover and/or black ball regions in the flight 

simulation) + 1; Input number of cloud cover or black ball 

regions required for the simulation 

Number of measured data points input 

NEX + 1 

NEX + 2 

Printer control character 

Parameter to eliminate the influence of the sun during the 

night portion of the flight simulation 

Parameter to eliminate the influence of the sun during the 

night portion of the flight simulation 

Array for storage of an 80 character title 

(Number of gas valvings during the flight +1) * 2; input number 

of valvings required for the flight simulation 

Number of data points collected for plotting 

NYC + 1 

NYC + 2 

Atmospheric pressure in subroutine ATM0S2 (millibars) 
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FORTRAN 

Variable 


Description 


PASL - Ambient air pressure at the launch site (mb) 

PATRO - Ambient air pressure at HTROP (mb) 

PATR02 - Ambient air pressure at HTR0P2 (mb) H 

PATR03 - Ambient air pressure at HTR0P3 (mb) j j 

PB - Balloon gas pressure (mb) j 

PB1(I) - Array of pressure altitude parameters corresponding to the array 
of times XPTT(I). PB1 » In (PASL/PB) 

PLOT - System subroutine used in plotting phase of program 

i 

PLOTS - System subroutine used in plotting phase of program 
PR - Prandtl Number 

j 

PTIMH - Time of day in hours ( Data are printed and saved for 

PTINC - Time increment in hours ( plotting in time increments of PTINC 

QF - Rate of heat transfer to the balloon film 

QG - Rate of heat transfer to the balloon gas 

R - Radius of balloon (Main program and subprograms CFH and CHRES) 

Universal gas constant (subprogram ATM0S2) 

REY - Reynold's number 
RGAS - Universal gas constant 

i 

RHOB - Air density at the base of the layers of the 1962 Standard ;; 

Altitude in subprogram ATM0S2 

RRES - Diameter'* of the balloon (subprogram CHRES) 

♦ 

RW - Infrared reflectivity of the balloon film 
RWSOL - Ultraviolet reflectivity of the balloon film 
SIGMA - Stefan-Bol tzmann constant 
SPHTG - Specific heat of the balloon gas 
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FORTRAN 

Variable 


Description 


SPHTW 

STOR 

TA 

TAFLO 

TAK 

TASL 

TATRO 

TATR02 - 

TATR03 

TAUW 

TB 

TBB 

TRBC(I) - 

TBDRP(I) - 
TCLDC(I) - 
TIME 


Specific heat of the balloon film 

Temporary storage location for the balloon gas viscosity 
Ambient air temperature (°K) 

Air temperature corresponding to altitude HFLOT in the four 
segment air temperature profile (°K) 

Ambient air temperature (°C) 

Air temperature at launch (°K) 

Air temperature corresponding to altitude HTROP in a four 
segment air temperature profile (°K) 

Air temperature corresponding to altitude HTR0P2 in a four 
segment air temperature profile (°K) 

Air temperature corresponding to altitude HTR0P3 in a four 
segment air temperature profile (°K) 

Infrared transmittance of the balloon film 
Air temperature at the base of the atmospheric layer of the 
1962 Standard Atmosphere in subprogram ATM0S2 
Black ball temperature (°K) 

Array of black ball temperatures corresponding to the cloud 
cover array ALBC(I) and the array of times TCLDC(I). Used 
when time exceeds TIMFL. (°K) 

Array of times associated with the ballast drops of the 
simulated flight (min.) 

Array of times associated with the cloud cover factors ALBC(I) 
and TBBC(I) (min.) 

Time of day at each point of the flight simulation (min) 
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FORTRAN 

Variable 


Description 


TIMFL 

TIMH 

TIMH1 

TIM 10 
TIMM 

TIMOL 

TIMSR 

TIMSS 

TOHIGH 


TR 

TRDEL 


TRK 

TS 

TSDEL 


Estimated time to float. Computed from ASCRT and HFLIT. 

Time of day at each point of the flight simulation (hrs) 

Time from launch at each point of the flight simulation 
(min) 

Temporary storage of time of launch (min) 

Time from launch of each point of the flight simulation 
(min) 

Time of launch (hrs. local standard time input; converted to 
minutes within program) 

Time of sunrise (min) 

Time of sunset (min) 

A parameter in subroutine ATM0S2 which indicates when the 
atmospheric pressure input is less than the capability 
of the subroutine and thus corresponds to an altitude 
that is too high to be accurately modeled by the sub- 
routine 

Temperature of the balloon gas (°K) 

Increment of the gas temperature during each step of the 
Runge-Kutta solution of the governing ordinary dif- 
ferential equations (°K) 

Temperature of the balloon gas (°C) 

Temperature of the balloon film (°K) 

Increment of the film temperature during each step of the 
Runge-Kutta solution of the governing ordinary differ- 
ential equations (°K) 
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FORTRAN 

Variable 

TSK 

TSL 

TT(I) 

TTIME 

TTMIN 

TTS 

TVALVE(I) 

TWILT 

TWSGl 

T1 

TIT 

T2T 

T3T 

T4T 


UMU 

VALVRT (I) - 


Description 

Temperature of the balloon film (°C) 

Maximum allowable change in balloon film temperature during 
each solution step 
See ETIME(I) 

Time from sunset or sunrise (min.) 

Plotting parameter. Minimum temperature to be plotted (°C) 

Time associated v'ith starting conditions (min.) 

An array of times corresponding to changes in valving rates 
(min) 

Length of twilight (hrs) 

Ultraviolet transmittance of the balloon film 
Temperature of the air in subprogram CFH (°K) 

Black ball temperature corresponding to altitude BALT in 
the black ball temperature profile (°K) 

Black ball temperature corresponding to altitude HBB2 in 
the black ball temperature profile (°K) 

Black ball temperature corresponding to altitude HBB3 in 
the black ball temperature profile (°K) 

Black ball temperature corresponding to altitude HFLOTin 

the black ball temperature profile (°K). Used for altitudes 
above HFLOT 

Viscosity of the balloon gas or ambient air, as appropriate 
Array of valving mass flow rates corresponding to the times 
TVALVE(I) (g/sec) 
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FORTRAN 

Variable 

VOLG 

VOLMAX 

X 


XLT 

XLT1 

XLT2 

XPTT(I) 


YPH(I) 
YPT1(I) 
YPT2( I) 

YPT3(I) 

ZSL 

Z1 


Description 

Volume of the balloon (m 3 ) 

Maximum volume of the balloon (m 3 ) 

Intermediate value in the computation of both the heat 
transfer coefficient between the film and gas 
(subprogram CHRES) and the heat transfer coefficient 
between the film and air (subprogram CFH) 

Physical length of the time axis on all plots (in) 

Physical spacing between sets of axes in the plotting (in) 
Physical spacing between sets of axes in the plotting (in) 
Array of times corresponding to the arrays of computed 
values of gas temperature, YPT3; air temperature, 

YPT1; film temperature, YPT2; balloon altitude, YPH; 
and the pressure altitude parameter, PB1 (min) 

Array of balloon altitudes, corresponding to the array 
of times XPTT(I), to be plotted(m) 

Array of air temperatures, corresponding to the array of 
times XPTT(I), to be plotted ( :J C) 

Array of balloon film temperatures, corresponding to the 
array of times XPTT(I), to be plotted (°C) 

Array of balloon gas temperatures, corresponding to the 
array of times XPTT(I), to be plotted (°C) 

Maximum allowable change in balloon altitude during each 
time step of the solution (m) 

Balloon altitude (m) 
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4 

Description 

21 DEL - Increment of balloon altitude during each time step of 

the solution (m) 

22 ** Balloon vertical velocity (m/min) l 

22DEL - Increment of balloon vertical velocity during each time 

step of the solution (m/min) ij 

i 

,i 


i’ 
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Table 5. Summary of Program Output 


Printed Output 

1. ) Heading 

2. ) Listing of input data in namelist FIPUT 

3. ) Tabular listing of the input ballast drop schedule: 

Time of ballast drop (sec) Ballast rate (g/min) 

(NBDRP values printed) 

4. ) Tabular listing of the input cloud cover details: 

Time of cloud cover (sec) Albedo Blackball Temp. (°K) 

(NCLDC values printed) 

5. ) Tabular listing of the input valving schedule: 

Time of the valving (sec) Valve rate (g/min) 

(NVALVE values printed) 

6. ) Tabular listing of the input measured data from the balloon 

flight being simulated: 

Time (sec. ) A1 titude(m) Air Temp.(°C) Film Temp; (°C) Gas Temp.(°C) 
(NEX values printed) 

7. ) The computed values of the effective emissivity of the balloon 

gas, the effective emissivity of the balloon film, and the effective 
interchange emissivity between the balloon gas and film. 

8. ) The computed values of the effective infrared absorptance of the 

balloon gas and the effective infrared absorptance of the balloon 
film. 
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ORIGINAL PAGE IS 

9.) Tabular listing of the computed results: POOR QUALITY 

Time(min ) , Altitude(m) , Velocity(m/min) , Air Temp. (°C) , 

Film Temp. (°C) , Balloon volume(m ) , Solar Flux , Blackball 
Temp. (°K) , Free lift (kg) . Pressure altitude (millibars) . 

Mass of gas (g) , Reynold's number . Coefficient of Drag , 

Payload mass (g) 

10.) The total number of data points saved for plotting (NYC). 

Plotted Output (Requires local subroutines equivalent to AXIS1, PLOT, PLOTS, LINE) 

1. ) A plot of the computed balloon altitude-time history and a plot of 

the measured balloon altitude-time history (if NEX is not 0), on a 
single set of axes. 

2. ) Plots of the computed balloon film and gas temperature and the ambient 

air temperature as functions of flight time. Also plots of measured values 
of film, gas and air temperature if NEX is not 0. 

3. ) A plot of the computed balloon pressure altitude parameter, 

ln < p launch site /P alt1tude ) » as a funct1on of time > and a similar P lot of 
the measured pressure altitude parameter if NEX is not 0. 
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Reference 


Carlson, L.A. and Horn, W.J., "A Unified Thermal and Vertical Trajectory 
Model for the Prediction of High Altitude Balloon Performance," Texas 
Engineering Experiment Station Report TAMRF-4217-81-02, June, 1981. 
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APPENDIX B 


SAMPLE CASE 

The input and output for a sample case is presented on the following 
pages. The 167N flight launched from Palestine, Texas on July 24, 1980, 
was chosen because it involves most of the features of the THERMTRAJ 
program with the exception of programmed valving. The plotted output 
should provide a good indication of the accuracy of the program in terms 
of predicted altitude as well as balloon gas and film temperature. 
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Table l.B. THERMtRAJ input data for the 167N flight 
(* Designates the column numbers of the data cards) 


, 33* 367050 1 21*567890 1 234 547090 12 3*56*890 123*3*7090 1 23496780 01 234 567 190 

* Fl1shT^ 167«N NUMBER 12 Si TH DOUBLE CD CN DESCENT AND HIGH CIRUS 

6FIPUT T[MaLalO.S8.CUTDT7.|8.S0,ASCnT-ZS7.3.t-RUlT*3583O.,PA5L-l‘}0«., 

M&SSFx i 94 2 E*5 • MG AS *ft • IE 4t MPAY I ■ I • 54fl2C5»VCL'MJ<»6 »6375C4» C AYL* 1 6 • t 

Al'l^i*l77l:fSiLl30A.55 HTROP- l 2500 . .TATRO-220.6 5.HTR0P2-1 3300 , . 

T A TR02*2 fll 35 «HTROP3*ZGOOO « » T ATRO 3*21 l .52 t HFUOT"35'300 • * T AFLO» 2 AS» 73* 

H 00 ?- 1 33 Q 0 • » H 803*2 C 909 . » T 2 T * 209 » A , T 37»209 .4 » TAT * 209 . A » DPLT * 1 • 0 « 
t ATM 5—6 * CEND 


0.1 0.0 

673 . 7257. *4 

676.125 0.0 

1200.0 ?»0 

0.0 A 0.2775 

170 . 0.2775 

5 * 5.0 0.13 

1 ? 00 . 3.18 
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204.4 

214.4 
214.4 
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Figure l.B. (continued) Altitude Trajectory for Flight 167N (66375 m 3 
balloon launched 1135 CDT, 24 July, 1980) 7 q 
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Figure l.B. (continued) Altitude Trajectory for Flight 167N (66375 
balloon launched 1135 CDT, 24 July 1980) 71 
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(continued) Altitude Trajectory for Flight 167N( G6375 m 3 
balloon launched 1135 CDT, 24 July, 1980) 72 
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Figure l.B. (continued) Altitude Trajectory for Flight 167N (66375 rrf 5 
balloon launched 1135 CDT, 24 July, 1980) 73 
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Figure 2.B. (continued) Temperature Profile for Flight 167N( 66375 m 
balloon launched 1135 CDT, 24 July, 1980) 
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Figure 2.B. (continued) Temperature Profile for Flight 167N (66375 m 
balloon launched 1135 CDT, 24 July, 1980) 
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Figure 2.B. (continued) Temperature Profile for Flight 

167N (66375 m 3 balloon launched 1135 CDT, 24 July 1980) 
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Figure 3.B. Pressure Trajectory for Flight 167N ( 66375 m 
balloon launched 1135 CDT, 24 July, 1980) 
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Figure 3.B. (continued) Pressure Trajectory for Flight 167N (66375 m 3 
balloon launched 1135 CDT, 24 July, 1900) 
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Figure 3.B. (continued) Pressure Trajectory for Flight 167N( 66375 m 
balloon launched 1135 CDT, 24 July, 1980) 
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Figure 3.B. (continued) Pressure Trajectory for Flight 167N(66375 m' 
balloon launched 1135 CDT, 24 July, 1980) 
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Figure 3.B. (continued) Pressure Trajectory for Flight 167N ( 66375 m 
balloon launched 1135 CDT, 24 July, 1980) 
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